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Abstract 


Flutter of an aeroelastic structure is potentially destructive aeroelastic instability. This phenomenon has 
motivated research within the aeroelastic community to develop methods that can accurately predict aeroelastic 
instabilities. The Flutterometer method used herein, and as developed by NASA DFRC, is based upon the // method 
which has been coupled with wavelet filtering processes in estimating aeroelastic models from flight data. The 
approach leads to a methodology to predict the occurrence of flutter boundaries, and may prove to reliably predict 
flutter boundaries during flight tests. An analytical model is used as the first estimate of the aeroelastic structural 
dynamics, and uncertainty operators are introduced into the system to model variations between the theoretical 
system and the physical system. The modelling uncertainties are then updated from experimental data. Although the 
model used did not work well with this particular experiment, a sensitivity analysis was additionally performed and 
improvements suggested. 


Introduction 


Aeroelasticity is the interaction of structural, inertial, and aerodynamic forces. Together, these forces may cause 
a structure, such as an aircraft wing, to become unstable and undergo aeroelastic instabilities such as flutter. Among 
the first to study this characteristic in detail was Theodorsen, and, along with Garrick*’^, he predicted flutter 
velocities and frequencies which were compared to experimental results. Several papers have been published on 
flutter, including methods for its suppression^ '^’^’^ The exact velocity at which flutter occurs depends on the physical 
properties of the structure and flow field around it. In the design of aircraft, the study of aeroelastic stability is 
required in order to determine a safe flight envelope. Many areas have been identified by the flutter community for 
developing an accurate flutter prediction method, including the time and cost associated with safely expanding the 
flight envelope in which no aeroelastic instabilities occur^’®. Several methods have been used in an attempt to predict 
flutter. Some of these methods are briefly discussed in Reference 9, 10, and 11 and are presented below. A 
traditional method analyses system parameters, such as damping to monitor aircraft stability. A prediction-error 
method was used to estimate damping errors, and extended to use Kalman filtering as well as time-domain and 
frequency-domain characteristics to estimate damping. But damping may be highly nonlinear as flight conditions 
vary; therefore, damping trends may be of limited usefulness, as stability may be indicated at a test point, but may 
not confidently be extrapolated to other flight conditions. Thus the flight envelope must be expanded slowly* 

Another method considers the interaction of two modes participating in the flutter mechanism to utilize a 
stability parameter that is varied quadratically with dynamic pressure'^ The technique considered several interacting 
modes and demonstrated a predictive method for higher order instabilities. It was limited though in applicability for 
general flight flutter testing because assuming only the coupling of a few modes and requiring their observation may 
be too restrictive 

Stability parameters have also been introduced in flutter margin methods that consider an autoregressive moving 
average process to describe aeroelastic dynamics While these margins are applicable to complex systems and 
require only turbulence for excitation, the flutter boundary is determined from the extrapolation of a nonlinear 
function and may be misleading. 


An approach currently being investigated is the use of a nominal model as the initial estimate of the properties of 
the aeroelastic structure and the incorporation of flight data in the model development process using uncertainty 
operators in a robust stability framework. The uncertainty operators are incorporated in the model through a linear 
fractional transformation (LFT), and the robust stability parameter known as the structured singular value^\ /i, is 
used to calculate flutter boundaries robust to these variations. In this manner, a worst-case flutter boundary is 
computed which accounts for the flight data directly. Difficulties of estimating a higher order model from data of 
low signal-to-noise ratio can be avoided, while time- varying dynamics are accounted for by updating the 
imcertainties between the nominal and true model The use of the structured singular value has led to this approach 
being dubbed “//-method”, which is the basis of a flutter-prediction method called the “Flutterometer”^’’’. The 
outcome may be immediately realized by on-line computing of robust flutter margins, with the result being the 
“distance” to the flutter boundary, determined by the largest increase in dynamic pressure (hence, the corresponding 
increase in velocity), for which the model is robustly stable. 


Background on li- Method Framework 

Any aeroelastic model is an approximation of the true system dynamics, but inaccuracies in various model 
parameters, as well as unmodelled dynamics must be considered in any stability analysis as well as control 
synthesis^ The characteristics of the system model are considered under the influence of perturbations, which 
contain the parameter uncertainties as well as any unmodelled dynamics or unmeasured forces^. The uncertainties 
are described by an operator, A , contained within a set, A , which is nonn bounded to reflect the limits of the 
perturbation sizes. Weighting matrices are often introduced in the // framework in order to normalize the uncertainty 
norm bound to unity. Thus, the perturbation set is defined as: 

A = {A:|HL<i} (1) 

By the Small Gain Theorem^’^\ the robust stability of a plant model P containing uncertainties is guaranteed for the 
set A if ||p||^ < 1 , though this condition is sufficient but not necessary because it may be overly conservative with 

respect to structured uncertainty. The equations of motion are rearranged into state-space form, and written as an 
upper LFT with the uncertainty block structure A . The LFT (Figure 1) is fonned with uncertainty over velocity, 



Figure 1. Linear Fractional Transformation block diagram with uncertainty matrix included. 


A^ , as well as parameter uncertainty, A^^^ , and dynamic uncertainty, , (which contains sensor and actuator 

uncertainties), thus achieving a robust aeroelastic modeP’^*’^^. The basis for the // method is the structured singular 
value, which is defined as: 
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and is a measure of the robustness of P with respect to A . The structured singular value, is a measure of 

the smallest destabilizing perturbation of P such that / — PA is singular for some A €. A. If no A 6 A exists that 
makes I — PA singular, then /y^(p):=0. In effect, the inverse of //^(P) is the magnitude of the smallest 

destabilizing perturbation A G A . For a unity norm bound set A , the system P is guaranteed to be robustly stable 
if (p) < 1 The structured singular value is computed numerically as an upper and lower bound using the /i-tools 

software of MATLAB®. 


Mode! Validation 


Robustness measurements are only meaningful if the uncertainty description of a given system model is realistic. 
Too much uncertainty may cause an overly conservative measure of robusmess, while too little uncertainty may 
generate a robustness measure that does not consider the true errors in the model^. It is necessary to validate any 
given model against a set of experimental data in order to determine a reasonable set of uncertainties. A /i-based 
method^’‘^’^‘ is used in this paper to validate the robust aeroelastic model of the airfoil pitch-plunge system, which 

assumes a known model, P , and its associated uncertainty set, A . An uncertainty set is sought such that (p, a) 

could generate a set of observed data, u (input) and y (output), in frequency space. The plant model is partitioned as 
follows: 
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where P^^ is the transfer ftmction from uncertainty input to uncertainty output, P,2 the actuator input to uncertainty 
output, P2, the uncertainty input to sensor output, and P^2 the actuator input to sensor output, which is the nominal 
plant transfer function. Define the following two matrices: 


^2 ” P \2^ 
hi ^Pii^-y 


(4) 


The model is not invalidated if the following is true: 

'4)>l (5) 

The above statement may sound counterintuitive, requiring the value of jj, to be greater than one to validate the data, 
while robust stability requires /i to be less than one. Consider the following relationship between u andy: 

0 = [P,,u -y] + P„A(7 - 

= 4+P„A(/-/>,A)-'P„ 

where the plant P - {Pwyhi^Pix^ hi } ' The existence of A satisfying (6) is equivalent to the plant not being 
robustly stable, thus, the model is not invalidated by the experimental data if ^(P) > 1 . In effect, this is an inverted 
robust stability criterion. 
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The validation algorithm makes use of three system matrices. The first is , which represents the true 
dynamics of the system and is used to validate a given set of experimental data to determine a realistic set of 
uncertainties. The second is , which is a theoretical approximation of the system dynamics (derived from the 
equations of motion), but does not take into account any unmeasured forces or unmodelled dynamics. This matrix is 
scaled with the uncertainties found using in the validation process to give the third matrix, P ^^^ , the system 

matrix which includes a perturbation in the freestream velocity and is utilized in the jd framework to determine the 
change in freestream velocity necessary to drive the system to instability. 


Equations of Motion for the 2 DOF Model 

The implementation of the Flutterometer is demonstrated using a wing section mounted to the Nonlinear 
Aeroelastic Test Apparatus (NATA)^^ in our wind tunnel (Figures 2 and 3). This system has been used extensively 
to study both linear and nonlinear aeroelastic motion, as well as develop control strategies'*’^'^'^'^’^^. The system 
consists of a wing section mounted on a carriage and is restricted in motion to pitch and plunge degrees of freedom. 
The bending and torsional modes of the wing are represented by springs and cams, which permit translational and 
rotational motion. A data acquisition system measures the pitch and plunge dynamics of the system, as well as 
provide control surface actuation, if applicable to the test. Figure 6 shows the schematic setup of the wing section 
with the data acquisition system. 
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Figure 2. The aeroelastic system permits two degree-of-freedom motion. A control surface is attached tothe 
trailing edge to provide controlled sinusoidal forces during tests. 



Figure 3. The model support system (NAT A) permits prescribed nonlinear motion in two degrees of freedom. 
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Figure 4. Data acquisition and active control schematic diagram for NATA tests. 


Flutter conditions may be reached on this system with minimal concern of structural failure, and parametric 
investigations of the system characteristics may be conducted, where the effect of changes on stiffness, elastic axis 
position, mass distributions, as well as various nonlinear features, may be explored. For this investigation, a 
linearized form of the equations of motion* for the wing section mounted to the pitch-plunge system is used: 


nyx^b 



(7) 


where h and a are the plunge and pitch coordinates respectively, nij^ is the total mass of the system, is the mass 
of the airfoil section, Xa is the nondimensionalized distance from the elastic axis to the center of mass, b is the semi- 
chord length, and are viscous damping terms, and are the pitch and plunge stiffness, and is the 
moment of inertia of the airfoil about the elastic axis. 


The aerodynamic lift and moment, L and M respectively, are given by: 
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where p is the air density, V is the freestream velocity of the wind tunnel, 5 is the airfoil span, a is the 
nondimensional distance from the midchord to the elastic axis, and are lift coefficients, and are 

moment coefficients, and P is the flap deflection. For this experiment, a frequency sweep is used for the control 
input. The values of all the parameters in the equations are listed below in Table 1 . 
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Table 1: System Parameters 


b 

0.135 m 


1 .664 kg 

s 

0.6 m 

nif 

0.273 kg 


2(0.3233)Z?-6-aZ^ 

m ^ 

0.114 kg 


r.Jb 

m„. 

m., +m„ 


2844.4 N/m 

nij. 

M + nif +m„ 

K 

3.5895 N.m/rad 

h 

0.00681 kg.m^ 

Ck 

27.43 N/m 

L 

0.001735 kg.m^ 

Ca 

0.036 kg.m^/s 

la 

0.04325f/,+/„ + «v'// 

P 

1 .225 kg/m^ 


3.358 

^ma 

6.28 
(0.5 + «)C,„ 


-0.635 


Here, , m^, and are the masses of the wing section without the control surface, the servomotor, and the control 
surface respectively, and and are the mass moments of inertia of the control surface and the servomotor 
respectively. 

The equations of motion are rearranged into state-space form and uncertainties in the parameters are introduced. 
A detailed derivation of the uncertain model can be found in Reference 26, introducing a perturbation in the 
freestream velocity, which is used to determine the flutter velocity through //-analysis, with the main points 
summarized here. The uncertain plant model corresponds to the following equation: 









= p 





_ y _ 


J ^ 


( 9 ) 


where , z^, z^, and z^ are input signals to the uncertainty block for freestream velocity, mass, damping and 

stiffness, and aerodynamic coefficient uncertainty, while and are the uncertainty feedback signals to 

the plant from the uncertainty block. Equation (9) can be written as an LET as shown in Figure 5. 
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Figure 5. LFT of the uncertainty model for the pitch-plunge system. 


where , A^, and A^ , are the blocks that contain the uncertainties in the parameters, 

is introduced to the freestream velocity, weighted by : 

The uncertainty, , 

u = u^+w^s^ 

(10) 

and is a convenient mathematical representation that allows operating conditions to be considered by setting S^j . 
Similarly, the uncertainties and weights for the mass variables are: 

»h =»ho 

m,. = ntur„ + 

and the damping and stiffness uncertainties: 

(11) 


Ca =^oo+KAa 

~ ^ho 

and finally, the aerodynamic parameters: 

(12) 



(13) 


Table 2 shows the initial uncertainty weight scaling factors used for each parameter uncertainty. An input 
multiplicative uncertainty has also been included in the uncertainty set. 
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Table 2; Uncertainty Weights for Associated Parameters 


ntj 

1 


1 


1 


2 


0.003 


140 


0.17 


1 


1 


1 


1 


0.2 


The Flutterometer 


The Flutterometer^, developed at NASA Dryden Flight Research Center using the MATLAB® /i-Tools software 
package, was used to determine flutter and instability margins for the wing section. The Flutterometer predicts flutter 
margins via a program that incorporates the /v-method to validate a system model against experimental data. Wavelet 
filtering processes have been incorporated with the /i-method for model identification.^^’^* Recently, the 
Flutterometer has been implemented on-line during flight tests on an F/A-18 research aircraft in which flight data at 
various test points was gathered and the flutter envelope determined^. The use of the F/A-18 explored the 
advantages of the Flutterometer over traditional methods®’^. Using the equations of motion that describe the response 
of the pitch-plunge system, was determined for the freestream velocity of each set of experimental data and an 
appropriate set of uncertainties was established. The model uncertainties were updated until the experimental data 
lay between the upper and lower bounds of the structured singular value. The second model, , was also 

determined as an approximation to the pitch-plunge system at a model velocity, U f^odei • Once an appropriate model 

was determined, P^^^ was determined from P^^^ by scaling P^^^ with the determined uncertainties and a variation 

was introduced to the freestream velocity as an additional uncertainty (as shown in Figure 5). Finally, the 
perturbation in the freestream was evaluated until the system was found to become unstable. 


Summary of Results 

Two methods were initially used to theoretically evaluate the flutter velocity: i) a classical V-g diagram 

approach, and ii) a robust stability analysis approach^ using the structured singular value, //. For both methods, an 
analytical model of the pitch-plunge system was used, determined from the equations of motion. For the first 
approach, an LFT was formed with the nominal model of the system at a given model velocity, U . An 

aeroelastic model was evaluated at different values of freestream velocities by varying (i.e., ~ 

where Sy is updated by a bisection search. At each new velocity, the eigenvalues of the system are evaluated until 

the velocity at which the system becomes unstable is reached. In turn, this velocity is used to determine the critical 
dynamic pressure at which flutter occurs. Figure 6 shows a typical plot of the linear pitch-plunge aeroelastic system, 
which would undergo flutter motion when reaching a destabilizing velocity. This plot was produced for 
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= 3.5895 N-m/rad and a = -0.68 . Beyond the critical velocity, the two aeroelastic frequencies coalesce thus 

indicating classical flutter that occurs as seen by the coupling of the two modes^. For different physical properties, 
another type of phenomenon is observed for the same system as shown in Figure 7, where a = -0.4 . Rather than 
under going flutter, the system undergoes aeroelastic divergence^^. The aeroelastic system fails, due to static 
instability, when it approaches this velocity. The freestream velocity at which divergence occurs is where one of the 
modal frequencies reaches zero. For approximately a < -0.6 , divergence is predicted; for a > -0.6 , classical 
flutter is predicted. We note that a = -0.5 represents the quarter-chord of the wing, which also represents the 
aerodynamic center and the center of pressure. 
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Figure 6. V-g plot for the pitch-plunge system indicating flutter for a - -0.68. 
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Figure 7. V-g plot for a pitch-plunge system approaching divergence velocity. 

The second approach makes use of an LFT of the system with the parametric and dynamic uncertainties, as well as 
. For a given model velocity, U , the uncertain system model is fonned. An initial Sy is chosen and the 

system is scaled by the uncertainties and by Wy5y , as described by equations (10)-(13). The structured singular 
value, /i, is evaluated for the scaled system and fVy is iteratively updated by fVy / /i to rescale the velocity 
perturbation (since the velocity is the perturbation parameter of interest, and, applying the definition of here, 
/j = l/(fVySy)). The process is repeated until /i<l. The predicted flutter velocities for each model are 

summarized in Table 3. For a = -0.68 and ~ 6 m/s, the V-g approach predicted a flutter velocity of 

15.1791 m/s while the jj, prediction gave 14.304 m/s. Note that for each model velocity, the predicted flutter 
velocities compare extremely well with the nominal case. No matter what model velocity was used initially, the 
same flutter velocity could be predicted. For any given initial model velocity, the predicted flutter velocities using 
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the V-g method agreed within the order of 0.0001 between models as shown in Table 3. While the same type of 
agreement is desired here, the results from the robust stability analysis differed a little more between models. This 
difference can be expected since the model is being approximated at each velocity, and uncertainties have been 
introduced. The flutter velocity predictions obtained from robust stability analysis are always less than the nominal 
predictions. 


Table 3: Predictions of Flutter Velocities for a=-0.68 for the Analytical Model 




U „ (m/s) 

6 

15.1791 

14.304 

7 

15.1793 

14.386 

8 

15.1796 

14.479 


Experiments were then performed for freestream velocities up to ~10 m/s, slowly approaching the flutter 
velocity, but not actually achieving it. Using simulated data, the Flutterometer predicted flutter velocities of 
approximately 13 m/s for the linear setup. The system itself has been observed to undergo flutter at a freestream 
velocity of ~14 m/s on the linear configuration. A command input to the trailing edge control surface was entered 
(with a steadily increasing frequency) and the response was measured. Initially, the command input was started at 1 
Hz and increased 0.5 Hz every 10 cycles for approximately 45 seconds. This frequency selection did not capture the 
modes between 1.6 and 2.3 Hz (where the natural frequencies of the pitch and plunge modes lie). The first 
frequency sweep time history is shown in the Figure 8. The frequency sweep was then modified to start at 1.5 Hz 
and increase by 0.2 Hz every 10 cycles. Figure 9 shows this time history. 
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Figure 8. Time history of frequency sweep. A 0.5 Hz change in frequency is provided every 10 cycles. 

The initial test frequency is 1 Hz. 

The experimental data was compared to a simulation from the linear model using the same input history. Some 
significant differences were noted between the two results, particularly in the pitch mode, as is shown in Figure 10, 
for a freestream velocity of 9.29 m/s. Interestingly, when the model response was multiplied by 2.5, both data sets 
compared nicely. It was not understood what caused the differences in amplitude, but the frequencies were in exact 
agreement. 
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Figure 9. Time history of frequency sweep. A 0.2 Hz change in frequency is provided every 10 cycles. 

The initial test frequency is 1.5 Hz. 


With the experimental data obtained, the Flutterometer showed no convergence (the fj, value initially computed 
was too high and would require a substantial number of iterations for the method to converge, if it w'as possible for 
this particular set of data, but was terminated early, with no results obtained). The amount of uncertainty required 
was too great to place the experimental data within the //-bounds. See Figures 1 1 and 12 for the time history of the 
pitch motion, with the model data included for comparison. Using this disagreeing experimental data in the 
Flutterometer resulted in a non-convergence of the algorithm to any flutter velocity, and an extremely high 
computational time. 
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Figure 10: Time response of pitch motion, comparing experimental to simulated response. 
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Figure 11: Time response of pitch motion, with model response multiplied by 2.5. 

A sensitivity analysis was performed to examine how much error can be introduced into various parameters to 
determine the quality of the flutter velocity predictions due to model errors. Uncertainties, one at a time, of 30 
percent were introduced into the stiffness, mass, viscous damping, and aerod)mamic coefficients, and simulated data 
was produced with this erroneous model. The Flutterometer predicted a flutter velocity of 1 1.7 m/s. In addition, a 
30 percent uncertainty was also introduced into all the aerodynamic coefficients at the same time. Again, the 
predicted flutter velocity was 1 1.7 m/s. The Flutterometer showed no sensitivity to the above parameter changes and 
converged to a flutter velocity. A second sensitivity test was also performed using simulated data by adding random 
noise at the input and the output to investigate how much error in the data the Flutterometer could tolerate. It was 
found that if the random noise level added to the simulated data was less than 0.0003, the Flutterometer converged in 
five iterations to a predicted flutter velocity of 11.7 m/s, else, no convergence was achieved. Note that the above 
tests were performed for a = -0.68 . 

One explanation for the anomalies seen in Figure 10 is the fact that the model used in the data validation process 
was valid for small angles of attack. The rather large pitch an^litudes may be affected by aerodynamic stall effects, 
which contribute to significant nonlinear responses not accounted for in the model. Stall effects are briefly described 
in Reference 30. This type of behaviour has also been observed in experiments described in Reference 5. 

A second explanation is that the commanded deflection for the control surface was not the actual deflection the 
control surface was undergoing. The wing section model used did not ha\'e any method for measuring the true 
control surface deflection, as the actual deflection would have been more appropriate to use as the input to the 
system. 


Concluding Remarks 

A flutter predicting methodology, the Flutterometer, was examined as a technique to predict the occurrence of 
flutter. Classical predictions were made using a V-g approach and a /^-analysis approach. By using the ffeestream 
velocity as an uncertainty parameter, a matched-point solution is achieved for the flutter velocity for the V-g 
approach, no matter what operating freestream velocity is selected for the model. Using the Flutterometer, an 
attempt was made to implement experimental data to the Flutterometer. An analytical model was determined for the 
wing section system at a particular velocity and the uncertainties updated. Unfortunately, since the experimental data 
did not agree well with the response produced by the model, particularly in the pitch degree of freedom, the result 
was a non-converging process, with high computational time. A sensitivity analysis showed that the Flutterometer 
algorithm worked well for small errors. Some of the observed responses, particularly the large pitch amplitudes, 
could not be adequately explained. In this case, improvements to the current model/experiment hardware associated 
with our two degree-of-ffeedom system should be examined. 
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